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The NJL model of dense three-flavor matter with axial anomaly: 
the low temperature critical point and BEC-BCS diquark crossover 
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We study the QCD phase structure in the three-flavor Nambu-Jona-Lasinio model, incorporat- 
ing the interplay between the chiral and diquark condensates induced by the axial anomaly. We 
demonstrate that for an appropriate range of parameters of the model, the interplay leads to the 
low temperature critical point in the phase structure predicted by a previous Ginzburg-Landau 
analysis. We also show that a Bose-Einstein condensate (BEC) of diquark molecules emerges in the 
intermediate density region, and as a result, a BEC-BCS crossover is realized with increasing quark 
chemical potential. 



PACS numbers: 12. 38. Aw. ll.10.Wx, 11.30.Rd, 03.75.Nt 



I. INTRODUCTION 

The phases of strongly interacting matter described by 
quantum chromodynamics (QCD) at finite temperature 
T and quark chemical potential (i is being actively stud- 
ied theoretically, as well as experimentally in ultrarela- 
tivistic heavy ion collisions at RHIC (Relativistic Heavy 
Ion Collider) and in the near future at the LHC (Large 
Hadron Collider). At low T and /i, the hadronic phase 
is realized with chiral symmetry dynamically broken by 
condensation of quark-antiquark pairs, the chiral con- 
densate (qq). On the other hand, at low T and high fi 
a color superconducting (CSC) phase characterized 
by formation of quark-quark pairs - a diquark conden- 
sate (qq) - is expected to appear owing to the attractive 
one-gluon exchange interaction or the instanton-induced 
interaction in the quark-quark channel. At high T for 
any fi, the quark-gluon plasma (QGP) phase |2J is real- 
ized with both the chiral and diquark condensates melted 
away. The phase transition from the hadronic phase to 
the QGP phase is indeed confirmed by recent lattice QCD 
Monte Carlo simulations indicating a smooth crossover 
for physical quark masses @, H| ■ 

Nevertheless, the first-principles lattice technique 
based on importance sampling is not applicable to QCD 
at finite \i due to the complex fermion determinant. 
This is why our understanding of the transition from 
the hadronic phase to the CSC phase relevant to the 
compact star physics is still immature and we have to 
basically rely on specific models of QCD, such as the 
Nambu-Jona-Lasinio (NJL) model [a, a|, the Polyakov- 
Nambu-Jona-Lasinio (PNJL) model |7H10l| , and the ran- 
dom matrix theory (RMT) [llH13j . These model studies 
together with the lattice QCD results have revealed the 
possible existence of the critical point [H], EH at high 
T between the hadronic phase and the QGP phase (see, 
however, [l6l|). 

Recently, we have pointed out the possibility of a new 
low temperature critical point between the hadron phase 
and the CSC phase in three-flavor QCD, on the ba- 
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FIG. 1. Schematic phase structure with two light (up and 
down) quarks and a medium heavy (strange) quark. In the 
hadronic phase, qq pairs condense, while in the color super- 
conducting (CSC) phase, the dominant condensation is qq 
pairing. In the quark-gluon plasma (QGP), all symmetries are 
restored without any pairing, while in the coexistence (COE) 
region qq and qq pairings coexist. The double line denotes a 
first-order phase transition. Adapted from Ref. [l8l |. 



sis of mo del- independent Ginzburg-Landau (GL) theory 
[ttI [l8j : the attraction between the chiral and diquark 
condensates induced by the axial anomaly leads to this 
critical point and an associated smooth crossover. Fig- 
ure [T] illustrates the dense three- flavor QCD phase di- 
agram with the new critical point at low T [18| . This 
may provide a mechanism of continuity between hadronic 
matter and quark matter (hadron-quark continuity) con- 
jectured by Schafer and Wilczek [l9j . Moreover, the 
idea of hadron-quark continuity is corroborated by recent 
studies on the spectral continuity of Nambu-Goldstone 
modes |18| and vector mesons [20(, and the formal sim- 
ilarity of the partition functions in the universal regime 
between the low and high fi regimes in three-flavor 
QCD at finite size, large compared with the inverse gap, 
but small compared with the pion Compton wavelength 
21]. In two-flavor QCD, similar new criti cal p oints have 
also been found in the NJL model |22h24| . although 
their origin is related to a repulsive vector-channel four- 
fermion interaction |22j , or electric charge neutrality and 
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/3-equilibrium conditions [23l . |24| rather than the axial 
anomaly. These studies may imply smooth crossovers 
not only as a function of T at low \x but also as a func- 
tion of n at low T in the realistic QCD phase diagram. 

The Ginzburg-Landau analysis for three- flavor QCD in 
Refs. [13, EH depends on the assumption that the mag- 
nitudes of the chiral and diquark condensates are suffi- 
ciently small near the phase boundaries, which may not 
be justified over the entire region in the QCD phase dia- 
gram and for strongly first-order phase transitions. The 
question is unanswered as to whether such a new critical 
point induced by the axial anomaly does really emerge 
in the (/z, T)-plane within the framework of phenomeno- 
logical models with reasonable parameters. 

The purposes of this paper are two-fold. First, using 
the three-flavor NJL model incorporating the interplay 
between the chiral and diquark condensates induced by 
the axial anomaly, we study the location of the new criti- 
cal point predicted in [17], [l8| . We demonstrate that this 
critical point indeed appears in the phase diagram for an 
appropriate range of parameters. Second, we show that 
the axial anomaly also triggers, in this model, a crossover 
between a Bose-Einstein condensed state (BEC) of di- 
quark pairing and Bardeen-Cooper-Schrieffer (BCS) di- 
quark pairing. This BEC-BCS crossover is identical in 
structure to that in nonrelativistic condensed matter sys- 
tems [25T - [27| ; and is discussed for relativistic systems [28l — 
[37j : the change in size of Cooper pairs at lower /x within 
the QCD Schwinger- Dyson approach (28l |. the relativis- 
tic BEC-BCS crossover of diquark pairing in the NJL- 
typ e model [29|, lU |33W36| and the diquark-quark model 
[30I [Hj], and a possible evolution from baryons in nu- 
clear matter to diquarks in quark matter with increasing 
/j, (37} are elucidated. Remarkably, as we will show in 
this paper, in relativistic quark matter at high /j, the ax- 
ial anomaly enhances the attractive interaction between 
quarks, leading to the emergence of a BEC state of di- 
quark pairing. 

To illustrate the essential physics induced by the axial 
anomaly and to avoid complications of charge neutral- 
ity and /3-equilibrium conditions, we assume SU(3) fla- 
vor symmetry m u = m<± = m s = m q throughout this 
paper. The generalization to include these effects will be 
reported elsewhere. 

This paper is organized as follows. In SecHU we formu- 
late the three-flavor NJL model incorporating the inter- 
play between the chiral and diquark condensates induced 
by the axial anomaly. In Sec. we discuss the phase 
structures with and without the interplay. In Sec. lIVl we 
show that the interplay leads not only to the new critical 
point but also to the BEC-BCS crossover of the diquark 
pairing at high density. Sec. |V]is devoted to a summary 
and concluding remarks. 



II. NJL MODEL WITH AXIAL ANOMALY 

The Lagrangian of the Nambu-Jona-Lasinio (NJL) 
model with three-flavors consists of three terms: 

C = g(* 7/1 0" -m q + fi lQ )q + £< 4 > + £( 6 \ (1) 

where q = (u, d, s) T is transpose of the quark field, m q is 
a flavor symmetric quark mass (m u = nrid = m s ). 
and are the four-fermion interaction and six-fermion 
interaction, respectively. The standard choice of is 

SI, 

£W=4 4 )+4 4) , (2) 

8 

4 4) =GE[(^) 2 + (^75T- a9 ) 2 ] 

a=0 

= 8Gtr(0V), (3) 

4 4) =# E [{q^T A \A'Cq T ){q T Ci lh T A \ A ,q) 
A , A '=2,5,7 

+ {qT A X A ,Cq T ){q T CT A \ A ,q)\ 
= 2Htr[d) L d L +d R d R ] 1 (4) 

where fcj = {qR) J a {qL) l a , {d L ) ai = e abc e l:jk (q L y b C{q L ) k , 
and (d R ) ai = tabc£iik{qR) 3 b C{q R ) k cl with a, b, c and k 
the color and flavor indices, and C the charge conjugation 
operator, tr is taken over the flavor indices. The flavor 
U(3) generators r a (a = 0, • • • ,8) are normalized so that 
tr[r a T{,] = 2<5 a fc, and r A and A^' with A, A' = 2,5,7 are 
antisymmetric generators of flavor and SU(3) color, re- 
spectively. The coupling constants G and H with dimen- 
sion (mass) -2 are assumed to be positive. Starting from 
the one-gluon exchange interaction and apply a simple 
Fierz transformation, we obtain the ratio H/G = 3/4. 
However, we treat G and H as independent parameters 
of the effective Lagrangian and as detailed in Sec. lIIII take 
the values common in the literature. 

The four-fermion interactions introduced above are in- 
variant under SU(3) l x SU(3) r x V(l) A x U(l) b symme- 
try. The interaction C^P produces attraction of qq pairs 
in the color-singlet and spin-parity channel, inducing 
dynamical breaking of chiral symmetry with formation of 

a chiral condensate (38|. Similarly leads to attrac- 
tion of qq pairs in the color- anti-triplet and spin-parity 
channel, inducing color-flavor locked (CFL) supercon- 
ductivity with formation of a diquark condensate j39j . 

The six-fermion interaction in our model consists of 
two parts, 

£( 6 )=4 6 >+4 6 j. (5) 

is the standard Kobayashi-Maskawa-'t Hooft (KMT) 
interaction (It! Ell] . 



4 6) = -8K (det <t> + h.c.) . 



(6) 
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This interaction, invariant under SU(3)l x SU(3)i{ x 
U(1)b symmetry but not under U(1),a symmetry, ac- 
counts for the axial anomaly in QCD due to instan- 
tons. For positive coupling constant K with dimen- 
sion (mass) -5 , as we assume, the rj' meson has a larger 
mass than the other pseudo Nambu-Goldstone bosons 
(ir, ij, K). The term © serves the role in the QCD phase 
structure of making the chiral phase transition first-order 
as a function of T at u = for massless three-flavor limit 

As pointed out in [13, LL8| > the instanton couples the di- 
quark condensate and the chiral condensate, which mod- 
ifies the QCD phase structure in the intermediate density 
region. The effective interaction between the chiral and 
diquark pairing fields is described by a six-fermion term, 



h.c. 



(7) 



which has SU(3)lxSU(3)_rxU(1)b symmetry but breaks 
U(1)a symmetry explicitly. It is this term that is respon- 
sible for the low temperature critical point. We assume 
K' > 0, so that qq pairs in the positive parity channel, 
(di) = —(da), are energetically favored, as suggested 
from the weak-coupling instanton calculations [43j, l44j . 
Since the term (J7J) acts as an external field for x, it washes 
out the first-order chiral phase transition at intermediate 
density for sufficiently large K'\{d R )\ 2 [T3, GJ]. If we 
start from the instanton vertex and apply a simple Fierz 
transformation, we obtain the ratio K' / K = 1. However, 
since there is no a priori reason that K and K' have this 
ratio in the effective Lagrangian level, we keep them as 
independent parameters. 

The favorable condensates by the interaction + 

£(6) 

are the flavor-symmetric chiral and diquark conden- 
sates in the spin-parity + channel, defined by 



X$ij = 
sS A A' = (q T Cj 5 T A X Al q}. 



(8) 
(9) 



Here the condensate order parameters x an d s ? which are 
proportional to the order parameters a and d defined in 
the previous Ginzburg-Landau analysis [l7l . are re- 
lated to the parameters <f> and 4fi defined below Eq. Q 
by 



xSij = 2(</>y), 

sSai — 2((d L )ai 



-2((d R )ai 



(10) 

(11) 



We work at the mean-field level, linearizing the the 
products of operators X and Y as X 2 — > 2(X)X — (X) 2 , 
XY -> (X)Y + (Y)X - (X)(Y), and X 2 Y -> (X) 2 Y + 
2(X)(Y)X - 2(X) 2 (Y). Subtraction of the constant 
terms avoids double counting the interactions. In mean- 
field deviations from factorization are partially compen- 
sated for by redefinition of the coupling constants G, H, 



K, and K' . Then £< 4 ) and reduce to 

£W -^AGxqq-eGx 2 , 

Cf -> H [s*(q T C l5 T A X A q) + h.c] - 3H\ 



4 6 ) ^~2K X 2 qq + 4K X 3 , 

T 



£ S "> -^l s | 2 99 ~ ?X [s*(q T C l5 r A X A q) + h. 
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3K \ ,2 

+ — |s| 2 X- 



(12) 



Here and below we implicitly sum over A = 2,5,7 unless 
otherwise stated. 

To derive the thermodynamic potential, it is most con- 
venient to work in the Nambu-Gor'kov formalism; we in- 
troduce the bispinor field 



(13) 



with q = Cq T (and q c — q T C) the charge-conjugate 
quark field. Then the linearized form of the NJL La- 
grangian becomes 



C = ^S 1 " 1 * - U. 



(14) 



Here S 1 (p) is the inverse propagator in the momentum 
space: 



M 



Aj 5 t a X a 



(15) 



v -A*j 5 t a \ a 7 M p^ - ^70 - M 
where the dynamical Dirac mass in the gq-channel reads 

M( X , s, m q ) = m q - 4 (g - ^K X ^j X + ~K'\s\ 2 ,(W) 

and the dynamical Majorana mass in the gg-channcl 
reads 



A(x,s) 



2 H 



1 



K' X 



(17) 



They are both dependent on the order parameters, x an d 
s. The constant term needed to subtract double counting 
of the interactions in C is 

U( X , s) = 6G X 2 + 3H\s\ 2 ~ 4X X 3 - ^K'\s\ 2 X . (18) 

The terms in Eqs. (|T6| . (ITTl) . and (H"8l) arc shown dia- 
grammatically in Figs. |2fa)-(d), (3Ja)-(b), and Sfa)-(d), 
respectively. The chiral-diquark coupling (the K'-tevm) 
enhances the attractions in both the qq and qq channels. 

The thermodynamic potential at temperature T and 
quark chemical potential /i is given by 



2tt 3 2 



+U(x,s),(19) 



where Tr is taken over the bispinor space with the factor 
1/2 in front to correct for double counting of degrees 
of freedom. Evaluating the trace and summing over the 
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FIG. 2. Four contributions to the Dirac mass M (constituent 
quark mass). The chiral condensate \ is denoted by a black 
circle, the diquark condensate s by a black square, and s* by 
a white square. 
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FIG. 4. Four contributions to the constant term U. 





(a) 



FIG. 3. Two contributions to the Majorana mass A (super- 
conducting gap). 



fermionic Matsubara frequencies p° — iu) n — (2n+l)iriT, 
we arrive at the thermodynamic potential [6|: 



n = 



where 



d 3 p 
(2~ 



2Tln(l + e^/ T )+ Wl ± l +%s), (20) 



(21) 
(22) 



are the dispersion relations for the quasiquarks in the 
octet and singlet representations, with E p = ^Jp 2 + M 2 , 



Ai = 2A, and A 8 = A. Equations (HU) and (J20J) imply 
that x < is energetically favored for non-zero m q . On 
the other hand, s is generally complex and the thermo- 
dynamic potential is a function of |s| 2 . 



III. PHASE STRUCTURE 

We now explore the effect of the attractive Z<T'-term in- 
duced by the axial anomaly on the phase structure in the 
(/i, T)-plane of the three-flavor NJL model. The phase 
structures can be determined numerically by looking for 
the values of x an d s that minimize the thermodynamic 
potential in Eq. (|20[) globally. We follow the parameter 
choice of [6( where the coupling constants G and K are 
chosen to fit empirical mesonic quantities and the chiral 
condensate in the QCD vacuum. Table U shows two sets 
of parameters we adopt below. We vary the strength of 
the chiral-diquark coupling (the K' term) by hand. In 
order to illustrate how the anomaly changes the conven- 
tional phase structure and to avoid the complications of 
charge neutrality and /3-equilibrium, we assume SU(3) 
flavor symmetry, m u = ma = m s = m q . 



A. Without the chiral-diquark interplay 

We first show the phase structures without the K'- 
term in Fig. [5] Panels (a) and (b) show the results of 
the case of massless quarks, I, and finite mass quarks, 
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m, [MeV] 


GA 2 


HA 2 


KA 5 


M [McV] 


X V3 [MeV] 


I 





1.926 


1.74 


12.36 


355.2 


-240.4 


II 


5.5 


1.918 


1.74 


12.36 


367.6 


-241.9 



TABLE I. Two sets of parameters in the present three-flavor 
NJL model: the current quark mass m q , coupling constants 
G, H, and K, with a spatial momentum cutoff A = 602.3 MeV 
6]. The dynamical quark mass M and the chiral condensate 
X in the vacuum are also given. 



II, respectively. The phase diagram contains of a normal 
(NOR) phase defined by x = s = 0, a Nambu-Goldstone 
(NG) phase defined by x and s = 0, and a color su- 
perconducting (CSC) phase defined by % = and s ^ 00 
The chiral phase transition between the NG and NOR 
(or NG and CSC) phases is first-order, with the chiral 
condensate x changing discontinuously, while the color 
superconducting phase transition between the CSC and 
NOR phases is second-order, with the diquark conden- 
sate s changing continuously but not smoothly with a 
discontinuity in the diquark susceptibility ds/dT. 

In case II, the current quark mass changes first-order 
chiral phase transition to a crossover at high tempera- 
ture, whereas the first-order transition at high density 
region still remains, as shown in Fig. [5jb) . As a result, 
the second-order critical point, the Asakawa-Yazaki point 
0, El, appears in the (fi, T)-plane. The QCD criti- 
cal point moves down towards the [i-axis with increasing 
quark mass m q . The region x ~ is characterized by 
explicit breaking of chiral symmetry by the quark mass, 
while in the region x 7^ 0, chiral symmetry is dynamically 
broken. 



B. With the chiral-diquark interplay 

When the strength of the chiral-diquark coupling due 
to axial anomaly, K' , is relatively small (K' < A.IK 
in case I and K' < 3.8K in case II), the topologies of 
the phase structures remain unchanged, as one sees in 
Figs.[5fa) and (b). On the other hand, once K' exceeds 
a critical value K' c , the topological structure of the phase 
diagram changes as seen in Fig.[5](shown for K' = 4.2K): 
as discussed in [l?], [H| using the Ginzburg-Landau ap- 
proach, the K'-term, which acts as an external field for 
X, turns the first-order chiral phase transition into a 
crossover, and leads to a low T critical point at interme- 
diate density. As a result, the coexistence (COE) phase 
defined by x and s ^ spreads over the higher den- 
sity region across the second-order phase boundary from 
the NG phase in both cases I and II. The emergence of 



Even when chiral symmetry is broken only slightly (\ ~ 0) by 
the current quark mass, we use the same classification in terms 
of NOR, NG and COE as in Fig. Ob). 
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(b) massive case 

FIG. 5. The phase structure in the {fx, T)-plane in the three- 
flavor NJL model without the axial anomaly for (a) case I, 
massless quarks, and (b) case II, finite mass quarks. Phase 
boundaries with a second-order transition are denoted by a 
single line and a first-order transition by a double line. The 
dashed-dot line at high T in panel (b) shows the effective 
chiral crossover line, at which the susceptibility dx/dT peaks. 
See the text for further detail. 



the COE phase is consistent with the model-independent 
result that the chiral condensate x is proportional to the 
instanton density (or the strength of the axial anomaly) 
in the CFL phase |4J|. 

In Fig. [71 we depict K' c as a function of K for several 
values of the current quark mass, m q — 0, m q = 5.5 
MeV, and m q = 140.7 MeV. The i^-line separates the 
crossover and first-order regions; the chiral-diquark cou- 
pling K' favors the crossover, while the triple chiral cou- 
pling K favors first-order. As m q increases, the crossover 
region is enlarged since the current quark mass acts as 
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(b) massive case 

FIG. 6. Phase structure in the (/j,, T)-plane in the three-flavor 
NJL model with the axial anomaly for (a) massless quarks, 
and (b) finite mass quarks. The phase boundaries with a 
second-order transition are denoted by a single line and a first- 
order transition by a double line. The BEC-BCS crossover 
(dotted) line in (a) and (b) is defined by /j, — M(fj,,T), the 
dynamical quark mass. 



an external field on the chiral condensate, weakening the 
chiral transition. 



IV. BEC-BCS CROSSOVER INDUCED BY THE 
AXIAL ANOMALY 

The axial anomaly, for sufficiently large chiral-diquark 
coupling K' , not only triggers the low T critical point, 
but also a BEC-BCS crossover in the COE phase, as 
discussed in [37j in qualitative analogy with the cold 
atomic gases in condensed-matter physics. Physically the 
BEC regime is characterized by quark-pair sizes small 
compared to the interparticle spacing, while in the BCS 




1 - - 

K' c line at T = 

J 1 1 

12 3 

K/K 

FIG. 7. Critical lines in the (K, if')-plane at T = for 
several values of the current quark mass m q (K A 5 = 12.36). 
Chiral phase transition is realized as a smooth crossover in the 
region above the corresponding line while it is of first-order 
below the line. 



regime the pair size is large compared with the interparti- 
cle spacing. The possibility of a BEC-BCS crossover in a 
color superconductor, and the presence of a BEC regime, 
was first pointed out in [28| by looking at the change 
in size of the pairs with density. As shown later within 
an NJL-type model such a BEC regime appears for suffi- 
ciently large pairing attraction, H, in the gq-channel (35|. 
The novel feature we stress here is that the axial anomaly 
helps to realize the BEC regime through its contribution 
to the effective qq coupling in (|I7|) . 

H'=H+±K'\ X \. (23) 

Although H ~ 0.9G alone is not enough to produce the 
diquark BEC (see Fig. [5]), the chiral-diquark coupling K' 
increases H' sufficiently for a BEC to develop (see Fig. [6J. 

Analytically, the distinction between the BEC and 
BCS regimes lies in the nature of the quasiparticle dis- 
persion relations, Eqs. ([21]) and (|%2"]l. For (j, > M, the 
minima of the dispersion relations are at nonzero mo- 
mentum p = \J [X 2 — M 2 , with excitation gaps Ai and 
Ag, a structure characteristic of the BCS regime. On 
the other hand, for \i < M, the minima of the disper- 
sion curves are at p = 0, a structure characteristic of the 
BEC regime (26[. Figures [5] (a) and (b) show the curve 
fi = M(/j,,T) as the dotted line in the COE region. A 
BEC of bound diquarks exists between the solid and dot- 
ted lines. (Note that at T = the dotted line ends at fi 
= 286.6 MeV in case I and at n= 297.8 MeV in case II, 
reflecting the decrease of M(p, T = 0) from its vacuum 
value, Table H) 

The structure of the crossover from BEC to BCS, at 
the NG-BEC and NOR-COE boundaries in Figure H is 



7 



most clearly defined in terms of the diquark correlation 
function 

Gd(t,x) = -4H 2 (T t [s a (t,x)s a (0,0)]), (24) 

with sa(t, x) = g T (r, x)C'75T y 4A J 4g(r, x) (no summation 
over A). In the random phase approximation (RPA), 
this correlation function in the complex frequency (z) 
plane at temperatures above the diquark condensation 
temperature T c is given by 



G^(z,q = 0) 



1 

Iff 7 



p<A 



d 3 p 

(2tt) 3 2(£ p T|u)T2 



(25) 

where /(e) = l/(e e ' T + 1) is the Fermi distribution func- 
tion and A is the ultraviolet cutoff. As we see from 
Eq. (|25p. Go has a branch cut on the real axis for 
z > 2{M — fj,) (as well as a branch cut for z < — 2(M + /i) 
from the antiparticle contribution). 

In the regime fi < M(/x,T), for sufficiently large H', 
Gb(z,0) has a pole on the real frequency axis for < 
z = M D (fj,,T) - 2fi < 2{M - with M D (fi,T) the 
mass of a bound diquark. The system undergoes a BEC 
condensation when, at a given temperature, this pole first 
reaches zero frequency, G^ 1 (0,0) = [U H, fig ; thus 
the condition 



2/i = M D (/x,T), 



(26) 



determines the NG-BEC boundary (Fig. [5]). With in- 
creasing /x the branch cut starting at 2(M — /i) even- 
tually reaches down to z = 0, at which point the pole 
at the origin begins to move to complex values in the 
second Riemann sheet. For fi > M(/j,,T), the condition 
GJ^ 1 (0,0) = defines the onset of BCS pairing, stud- 
ied in detail in Ref. [46], and determines the NOR-COE 
boundary (Fig. [5]). 

Due to the abrupt change of x across the first-order 
transition line, the NG-BEC and NOR-COE boundaries 
are not smoothly connected; the former touches the first- 
order line at higher temperature Jhan the latter (see 
Fig.[B]). As discussed more fully in 18], this difference can 
be understood by noting that a larger chiral condensate 
X reduces the density of states at the Fermi surface and 
simultaneously increases the effective qq coupling ff ' [see 
Eq. (1231) ]. Since the latter effect dominates in the present 
parameter set, the critical temperature for diquark pair- 
ing is larger on the left side of the double line. 

The behavior of the Dirac mass M, the Majorana mass 
or gap A, and the bound diquark mass Mjj associated 
with the phase diagram in Fig.^b), are plotted in Fig. [5] 
as functions of [i for three temperatures, T=0, 62, and 80 
MeV. At zero temperature, we see the successive transi- 
tions from the NG phase, the BEC regime in the COE 
phase, to the BCS regime in the COE phase, with increas- 
ing fj,. The onset of BEC at fi = 279.2 MeV is determined 
by the condition M D - 2[i = 0. For T — 62 MeV, we see 
rather the successive transitions from the NG phase, the 
BEC regime in the COE phase, the NOR phase, to the 
BCS regime in the COE phase, with increasing /i. The 
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FIG. 8. The Dirac mass (M) and the Majorana mass (A) as 
functions of fi for T = MeV (top), 62 MeV (middle), and 
80 MeV (bottom). The excitation gap of the bound diquark 
in the medium, Mn — 2/x, extracted from the isolated zero of 
Eq. (|25|) . is also shown. The parameters are the same as in 
Fig. Ob). 



transition from the BEC regime in the COE phase to the 
NOR phase is first-order at fi = 284.5 MeV, with both M 
and A jumping discontinuously. The two phase transi- 
tions, from the NG phase to the BEC regime in the COE 
phase, and from the NOR phase to the BCS regime in 
the COE phase, signal the onset of a non-zero Majorana 
mass. Both transitions are correctly described in terms 
of the diquark correlation function [see Eq. ([2"5"|) ]. For 
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FIG. 9. The phase diagram in the (/i, 7f')-plane at T = for 
massless quarks, with the NG and COE phases. The BEC- 
BCS crossover in the COE phase for large K' is shown as a 
dotted line. The critical point and the critical end point are 
denoted by P and Q, respectively. 



T = 80 MeV, the system undergoes a first-order transi- 
tion, from the NG phase to the NOR phase at [i = 260.7 
MeV, which takes place before bound diquarks start to 
condense. 

Finally we consider the effect of the chiral-diquark cou- 
pling K' on the phases at T = 0. Figure [5] shows the 
phase diagram in the (/i, i"C')-plane for massless quarks. 
(A similar structure holds for finite mass quarks.) For 
small K', the system has an NG phase and a COE phase 
separated by a first-order line indicated by the double 
line which eventually terminates for large K' at the crit- 
ical point P. On the other hand, for K' sufficiently large 
compared with the cubic coupling, K, of the chiral field, a 
BEC regime of bound diquarks appears across a second- 
order phase transition (solid line) from the NG phase at a 
critical chemical potential /i = Mjj/2; the phase bound- 
ary joins the first-order line at the critical end point 
Q. The dotted line, /i = M(jjl,T), shows the BEC-BCS 
crossover; for somewhat smaller K' ', a novel first-order 
transition from the BEC to BCS regimes appears be- 
tween P and Q, with discontinuous changes of both the 
chiral and diquark condensates. 

V. DISCUSSION 

We have explored here the phase structure of dense 
three-flavor matter using the Nambu-Jona-Lasinio model 



incorporating the attraction between the chiral and di- 
quark condensates induced by the axial anomaly. We 
demonstrated that the low temperature critical point be- 
tween the hadronic phase and the color superconduct- 
ing phase predicted by the previous Ginzburg-Landau 
analysis [13, [H| indeed appears in the phase diagram for 
sufficiently large chiral-diquark coupling. We have also 
shown in Eq. (|23[) that the axial anomaly enhances the at- 
tractive interaction between quarks, leading to the emer- 
gence of a Bose-Einstein condensate (BEC) of diquark 
molecules. As a result, a BEC-BCS crossover in the di- 
quark pairing appears in the coexistence phase, which 
has both non-zero chiral and diquark condensates. 

In the phase diagram of the NJL model derived here, 
the BEC regime is realized adjacent to the lower den- 
sity Nambu-Goldstone phase of massive quarks. In QCD, 
however, the low density phase is in reality nuclear mat- 
ter; There remains the important problem of learning 
how the gas of bound diquarks and unpaired quarks un- 
dergoes a transition to a gas of three-quark bound states, 
or nucleons, at low density. Describing this transition will 
require going beyond the mean-field approximation, and 
RPA, to take into account residual interactions between 
the diquarks and unpaired quarks. Recent work [47l [48[ 
on mixtures of bosonic and fermionic atoms indicates a 
phase diagram very reminiscent of this scenario in QCD. 

It is also important to make our phase diagram more 
realistic by including effects such as the Fermi momen- 
tum mismatch induced by a strange quark mass, charge 
neutrality and /3-equilibrium. Open questions include 
whether the low temperature critical point can survive 
in an inhomogeneous chiral crystalline phase 

MM 

or 

the Fulde-Ferrell-Larkin-Ovchinnikov phase |5ll456l |. and 
how the phase structure obtained here is affected by 
the confinement-deconfinement phase transition charac- 
terized by the Polyakov loop [7H10| . We defer these prob- 
lems to future publications. 
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